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Surfaces and interfaces offer new possibilities 
for tailoring the many-body interactions that 
dominate the electrical and thermal properties 
of transition metal oxides mm- Here, we use 
the prototypical two-dimensional electron liquid 
(2DEL) at the SrTiO3(001) surface [4HZ1 to re¬ 
veal a remarkably complex evolution of electron- 
phonon coupling with the tunable carrier den¬ 
sity of this system. At low density, where super¬ 
conductivity is found in the analogous 2DEL at 
the LaA 103 /SrTi 03 interface [8l- [13| . our angle- 
resolved photoemission data show replica bands 
separated by 100 meV from the main bands. This 
is a hallmark of a coherent polaronic liquid and 
implies long-range coupling to a single longitu¬ 
dinal optical phonon branch. In the overdoped 
regime the preferential coupling to this branch 
decreases and the 2DEL undergoes a crossover to 
a more conventional metallic state with weaker 
short-range electron-phonon interaction. These 
results place constraints on the theoretical de¬ 
scription of superconductivity and allow for a uni¬ 
fied understanding of the transport properties in 
SrTi 03 -based 2DELs. 

Carrier concentration is a key parameter defining the 
ground state of correlated electron systems. At the 
LaA103/SrTi03 interface, the 2 DEL density can be tai¬ 
lored by field-effect gating. As the system is depleted 
of carriers, its ground state evolves from a high mobility 
2 DEL [ 4 ], into a two-dimensional superconductor iHini 
with pseudogap behaviour m and possible pairing above 
Tc[l 2 ]. An analogous 2 DEL can be induced by doping 
the ( 001 ) surface of SrTi03. As for the interface, the 


surface 2 DEL is confined by a band bending potential in 
SrTi03 and consists of an orbitally polarized ladder of 
quantum confined Ti t2g electrons that are highly mo¬ 
bile in the surface plane HdHH. Thus far, the surface 
2 DEL has only been studied at carrier densities around 
2 X 10 ^^ cm“^, approximately a factor of five higher than 
typically observed at the LaA103/SrTi03 interface [ 5 ]-[ 7 ]. 
In the following, we present new ARPES data extending 
to lower carrier densities that are directly comparable to 
the LaA103/SrTi03 interface. We achieve this by prepar¬ 
ing SrTiO3(001) wafers in situ, which results in well or¬ 
dered clean surfaces that can be studied by ARPES over 
extended timescales as they are less susceptible to the 
UV induced formation of charged oxygen vacancies re¬ 
ported for cleaved SrTi03 0 [3 [m [IS]. Details of the 
sample preparation are given in the methods section. 

Fig. [^a shows an energy-momentum intensity map for 
a 2 DEL with a carrier density of 7220 ~ 2.9 - 10 ^^ cm“^ es¬ 
timated from the Luttinger volume of the first light sub¬ 
band and the two equivalent heavy subbands (see Sup¬ 
plementary Information section 2 ). The most striking 
feature of this data are replica bands at higher binding 
energy following the dispersion of the primary quasipar¬ 
ticle (QP) bands. The replica bands are all separated by 

100 meV and progressively loose intensity but can be 
visualized up to the third replica in the curvature plot 
shown in Eig. b. Erom the equal energy spacing of 
the replica bands we can rule out that they represent 
distinct subbands arising from quantum confinement 
13 US]. We thus interpret the replicas as shake-off ex¬ 
citations involving a single non-dispersive bosonic mode 
coupling electrons of momentum k and k + q. Erom 
its energy of ~ 100 meV, we can identify this mode as 
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FIG. 1. A two-dimensional liquid of large polarons in SrTiOa. a,b Energy-momentum intensity map and curvature 
plot for a 2DEL with n 2 D — 2.9 • 10^^ cm“^ taken at a photon energy of 44 eV with s-polarization. Note the dispersive 
replica bands at higher binding energy arising from strong coupling to the LO 4 phonon branch of SrTiOs- c,d Data taken on 
a sample with n 2 D — 3.7 • 10^^ cm“^ using s and p polarized light with hiy = 85 eV to selectively excite the light xy and heavy 
xzjyz orbitals, e Dispersion of the main and first replica bands extracted from the data in c,d. f Photon energy dependent 
measurements showing the lack of dispersion along kz. All data were measured in the second Brillouin zone to avoid the 
minimum of the matrix elements at normal emission. 


the highest frequency longitudinal optical phonon branch 
(LO4) of SrTiOs [n]. A plasmon mode in the same en¬ 
ergy range HZIIIH] can be excluded from the negligible 
density dependence of the mode frequency observed in 
Fig.H Our ARPES data show that the coupling to this 
mode largely preserves the dispersion in the replica bands 
and thus must be restricted to small values of q. This 
is a hallmark of Frohlich polarons, quasiparticles formed 
by an excess electron dressed by a polarization cloud ex¬ 
tending over several lattice sites that follows the charge 
as it propagates through the crystal [HIIMI]. Such a 
large polaron state preserves band-like transport but has 
an increased effective mass m*. From our ARPES data 
we can directly quantify this effect. Using a parabolic fit 
of the band dispersion in Fig.j^a we find ^ 1.4 me- 
Using a bare band mass of mo = 0.6 mg mm this corre¬ 
sponds to a mass enhancement m*y/mQ ~ 2.3, indicative 
of an intermediate e-p coupling strength. 

In Fig. c-f we investigate a sample with slightly 
higher carrier density. Using s and p polarized light, 
respectively, we selectively excite electrons from the light 
xy band with ^ 50 meV occupied bandwidth and a shal¬ 
lower heavy band derived from xz/yz orbitals. Both 
bands show a clear peak-dip-hump line shape with a dis¬ 
persive replica band as summarized in Fig. e. This 
implies that the entire orbitally polarized 2DEL includ¬ 
ing the heavy xz/yz states, which are believed to be im¬ 
portant for superconductivity mm, is a polaronic liq¬ 
uid at low carrier densities. The lifting of the orbital 
degeneracy by approximately 20 meV can be attributed 


to quantum confinement in the band bending potential 
which increases the energy of out-of-plane orbitals [71 122] . 
The reduced dimensionality arising from quantum con¬ 
finement is visualized directly in Fig. [^f where we show 
the absence of dispersion along the surface normal kz, 
characteristic of two-dimensional electronic states. 

We now turn our attention to the strength of e-p cou¬ 
pling. Its systematic evolution with carrier density in 
the SrTiOs (001) 2DEL is shown in Fig. Using s- 
polarized light, we resolve a single dxy QP band at low 
density. Energy distribution curves extracted at kp show 
a strongly reduced weight of the coherent QP and at least 
two phonon satellites each separated by ~ 100 meV. 
Empirically, we find that the experimental spectra are 
well described by a Franck-Condon model using a single 
phonon mode of ^ 100 meV. Highly restricted fits using 
the characteristic Poisson distribution I^/Iqp = a//' jn\ 
for the intensity ratio of the n-th phonon satellite and 
the QP peak in the Franck-Condon model are shown in 
Fig.[^g-1. Details of the fits are described in Supplemen¬ 
tary Information. From this analysis we infer a quasipar¬ 
ticle residue Z 0.2 at the lowest density studied here. 
This is beyond the validity of the perturbation theory 
result for the Frohlich model of Z = 1 — (a/2 [23l 1^ . 
where a is the dimensionless coupling constant, placing 
the SrTiOs 2DEL in a theoretically challenging regime 
of intermediate coupling. Moreover, the ratio of lattice 
energy to electron kinetic energy is neither small nor in¬ 
dependent of doping, which excludes a Migdal-Eliashberg 
approach. We therefore estimate a from the experimen- 
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FIG. 2. Evolution of the 2DEL spectral function with carrier concentration, a-f, Raw energy-momentum intensity 
maps of 2DELs with increasing carrier concentration indicated in units of cm“^. g-1, Energy distribution curves at the Eermi 
wave vector indicated by a dashed white line in the corresponding image plots. An exponential background describing the tail 
of the in-gap state at ~ —1.1 eV has been subtracted from the raw EDCs (see Supplementary Information). In g-j we show 
hts to a Eranck-Condon model with a single phonon mode. The coherent quasiparticle and incoherent phonon satellites are 
coloured in yellow and blue, respectively. 1 shows a calculation of the spectral function using the conventional Eliashberg e-p 
self-energy for the high-density limit found in Ref. [7] and reproduced in the inset. 


tally determined quasiparticle residue Z using the results 
of diagrammatic quantum Monte Carlo simulations of the 
Frohlich model reported in Ref. [24]. This gives a ~ 2.8 
at the lowest carrier density, comparable to a = 2 — 3 re¬ 
ported for lightly doped bulk SrTiOs based on an analysis 
of optical conductivity data [25l [26| . 

As the carrier density increases, the Fermi wave vec¬ 
tor increases monotonically, and new subbands become 
discernible. Concomitantly, the spectral weight of the 
phonon satellites weakens and for densities above n 2 ~ 
9 • 10^^ cm“^ they can no longer be resolved experimen¬ 
tally. At high density e-p coupling is not only weaker but 
also of fundamentally different nature. This is illustrated 
in Fig. [^f,L In this regime, the quasiparticle dispersion 
shows a weak kink at an energy of ~ 30 meV and no 
signs of replica bands can be discerned in the raw data 
or in curvature plots (see Supplementary Information, 
Fig. S4), providing direct evidence for a suppression of 
the long-range Frohlich interaction. Instead, consistent 
with Ref. [ 7 ], we find that the spectral function and e- 
p self-energy of the high-density 2DEL can be described 
by Migdal-Eliashberg theory with A ^ 0.7 and the same 
coupling to the entire phonon density of states. 

In order to track the crossover from the polaronic state 
to short-range e-p coupling, we analyse the quasiparticle 
residue and effective mass of the dispersive bands close 
to the Fermi level as a function of 2D carrier concen¬ 
tration. The results are summarized in Fig. [^ where 
we plot m* obtained from fits to the QP dispersion and 


Zikp) from fits to a Franck-Condon model as described 
above. Three distinct regimes can be identified. Below 
ni ^ 4-10^^ cm“^ both quantities depend weakly on car¬ 
rier density and appear to saturate around m* 1.4 mg 
and Zikp) ~ 0 . 2 , respectively. At intermediate carrier 
concentrations of n 2 D = 4 — 9 • 10^^ cm“^, the polaronic 
state persists, as is seen most directly from the spectra 
in Fig. h-j showing a clear phonon satellite over this 
entire regime. However, approaching n 2 ~ 9 • 10^^ cm“^, 
where the phonon satellites are no longer resolved ex¬ 
perimentally, the quasiparticle residue increases by more 
than a factor of two to Z « 0.5 and the effective cou¬ 
pling strength decreases to a ~ 1.3. This is opposite 
to the trend expected for short-range e-p coupling m 
and thus strongly supports our identification of long- 
range interactions as described by the Frohlich model. 
We note that the weak coupling to the LO 4 branch close 
to 712 reported in Fig. has recently been confirmed 
for the LaAlOs/SrTiOs interface 2 DEL in a soft x-ray 
ARPES study reporting Z ^ 0.4 for a sample with 
7120 ~ 8 • 10 ^^ cm“^, in excellent agreement with our 
findings [28] . 

The effective mass decreases more slowly than the 
quasiparticle residue and saturates at m* « 1.0 mg at 
high density. This slow decrease is characteristic of the 
Frohlich model where Z < mo/m*. For weak to interme¬ 
diate coupling, the effective mass of Froehlich polarons 
can be approximated as m*/mo = 1/(1 — (a/ 6 ) [191 l2l|- 
As shown in Supplementary Information Fig. S 6 , this re- 
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FIG. 3. Effective mass and quasiparticle residue in the 
SrTiOa 2DEL. Evolution of the effective mass m* (blue sym¬ 
bols) and quasiparticle residue Z (red symbols) with carrier 
density. Different symbols indicate data taken on substrates 
annealed at different temperature. Closed red symbols are 
obtained from Franck-Condon fits, while the last value with 
open symbol in the adiabatic Migdal-Eliashberg regime has 
been calculated from Z — mQlm*. Error bars indicate the 
reproducibility of our results. An additional systematic error 
cannot be excluded. The background color encodes the bare 
band width of the first light subband calculated from the ex¬ 
perimentally determined kp shown in the top-axis, assuming 
a bare mass of ttiq = 0.6 mg. Dashed lines are guides to the 
eye. The dome shaped superconducting phase observed at the 
LaAlOa/SrTiOa interface is indicated in grey. 


lation systematically underestimates the effective masses 
obtained directly from the quasiparticle dispersion but 
reproduces their trend as a function of density. We ten¬ 
tatively assign this behavior to the effect of electron- 
electron interactions, which is not fully included in the 
analysis of Z in the polaronic regime. Indeed, the fac¬ 
tor between the two effective masses is approximately 
1.3, which is consistent with the mass-enhancement due 
to electronic correlations estimated in Ref. [7]. For den¬ 
sities above n 2 , described by Migdal-Eliashberg theory, 
Z = m/m* suggesting that the quasiparticle residue sat¬ 
urates around Z = 0.6 in the high-density 2DEL. 

We attribute the breakdown of the polaronic state 
at high carrier density to improved electronic screen¬ 
ing suppressing the long range Frohlich interaction. In 
order to estimate the cross-over from dielectric screen¬ 
ing in the polaronic regime to predominantly elec¬ 
tronic screening, we compare the polaron radius = 
(fi/2m*flLO,4)~^^^ mi with the electronic screening 
length. Using the experimentally determined parameters 
we find ^ 6 A. Given the relatively large dielectric 
constant of doped SrTiOs, electronic screening is in the 
Thomas-Fermi regime and the screening length can be 
estimated from ttf = ^ elec¬ 

tron liquid. Assuming a uniform charge distribution over 


a 2DEL thickness of 3 unit cells and a static dielectric 
constant eo = 100 , which provides a good description 
of the band-bending potential in the doped surface re¬ 
gion mm, we estimate that ttf = for a density 
n 2 D ~ 5 • 10^^ cm“^, slightly above ni where Z starts to 
increase progressively resulting in a reduced effective cou¬ 
pling constant a. At higher carrier densities, electronic 
screening will rapidly become more important as eo de¬ 
creases simultaneously with the increase in n. We note 
that the moderate electronic screening length ttf justi¬ 
fies the above use of an expression for 3D electron liquids. 
This basic picture is confirmed by our calculations of the 
spectral function reported in supplementary Fig. S5. Us¬ 
ing exact diagonalization of a model Hamiltonian with 
constant coupling parameter, we find pronounced replica 
bands at low carrier density and a rapid suppression of 
the replicas as the carrier density is increased over the 
relevant range. These calculations thus reproduce the 
key experimental findings of our study and support the 
idea that increasing electronic screening drives the ob¬ 
served breakdown of the polaronic liquid at high carrier 
densities. 

The above results demonstrate that e-p interaction in 
SrTiOs based 2 DELs is remarkably complex and strongly 
dependent on carrier density. This provides new in¬ 
sight into the superconducting pairing mechanism at the 
LaAlOs/SrTiOs interface, which has so far eluded ex¬ 
perimental investigation. As shown in Fig. supercon¬ 
ductivity in LaAlOs/SrTiOs interface 2 DELs [S UTOl [29] 
occurs in the polaronic low-density regime and its sup¬ 
pression on the overdoped side coincides with a gradu¬ 
ally decreasing coupling to the LO 4 branch. This sup¬ 
ports the notion that superconductivity is phonon medi¬ 
ated [13 EDI EH and suggests that the pairing potential 
is dominated by the exchange of high-frequency longi¬ 
tudinal phonons rather than soft modes. We note that 
our experimental spectral functions do not exclude the 
formation of large bipolarons, which have been discussed 
early on in the context of superconductivity in polar ox¬ 
ides [32l |33] . Superconducting susceptibilities calculated 
within our exact diagonalization scheme reported in sup¬ 
plementary information section 5 provide additional in¬ 
sight. Although limitations in the Hilbert space size pro¬ 
hibit a quantitative comparison with experiment, the cal¬ 
culations clearly show that the dominant pairing chan¬ 
nel has 5- wave symmetry. Moreover, the saturation of 
the superconducting susceptibility with carrier density 
found in this model provides evidence for a competition 
between the opposite trends of density of states and effec¬ 
tive e-p coupling underlying dome-shaped superconduc¬ 
tivity in SrTiOs 2 DELs. We note that coupling to the 
LO 4 branch of SrTiOs with comparable strength was also 
invoked to explain the anomalously high superconduct¬ 
ing critical temperature of FeSe monolayers on SrTiOs 
substrates laiai]. 

Methods The Nb-doped (0.5 wt%) and La-doped 
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(0.075 wt%) SrTiO3(001) surfaces were prepared by mild 
Ar+ sputtering (600 eV, 2 /iA, 5 min) followed by an¬ 
nealing in 2 • 10“^ mbar oxygen for 0.5 h at temperatures 
varying from 700 to 1000° C as monitored by an infrared 
pyrometer. Consistent with measurements on in situ pre¬ 
pared and cleaved SrTiOs, we find no spectral weight in 
the entire bulk band gap on freshly prepared surfaces (see 
Supplementary Information, Fig. 1) [5l[T5j[T6]. Quantum 
confined metallic states appear concomitant with a local¬ 
ized in-gap state at higher energy after exposure of the 
surface to synchrotron light O [in Hn . We find that the 
rate at which the 2DEL carrier density increases under 
the synchrotron beam increases with the annealing tem¬ 
perature and exploit this to stabilize low-carrier densities 
over the extended exposure times required for detailed 
ARPES measurements. ARPES measurements were per¬ 
formed at the SIS beamline of the Swiss Light Source, the 
105 beamline of Diamond Light Source and the 1^ beam- 
line of BESSY II at the Helmholtz-Zentrum Berlin. Data 
was acquired at T « 20 K with hu = 30 — 100 eV and en¬ 
ergy and angular resolutions of 10 — 25 meV and ~ 0.2°, 
respectively. 
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